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ABSTRACT 



Context. Due to observational constraints, dark matter determinations in nearby clusters based on weak lensing are still extremely 
rare, in spite of their importance for the determination of cluster properties independent of other methods. 
Aims. We present a weak lensing study of the Coma cluster (redshift~ 0.024) based on deep images obtained at the CFHT 
Methods. After obtaining photometric redshifts for the galaxies in our field based on deep images in the u' (1x1 deg"), and in the B, 
V, R and I bands (42'x52'), allowing us to eliminate foreground galaxies, we apply weak lensing calculations on shape measurements 
performed in the u* image. 

Results. We derive a map of the mass distribution in Coma, as well as the radial shear profile, and the mass and concentration 
parameter at various radii. We obtain M200C = 5.1^2 iX 10'"* Iij^Mq and C2ooc = 5.0;^2 5' good agreement with previous measurements. 
Conclusions. With deep wide field images it is now possible to analyze nearby clusters with weak lensing techniques, thus opening a 
broad new field of investigation. 
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' 1. Introduction 

Although weak lensing studies are of great importance in cos- 
mology, in particular concerning various aspects of the analy- 

' sis of large scale structures in the universe (e.g. Mellier 1999; 
Bartelmann & Schneider 2001 ; Schneider 2006), dark matter de- 
terminations in clusters based on weak lensing have been lim- 
ited to relatively distant objects, usually with redshifts larger 
than about 0.15. This is mainly due to the fact that the clus- 
ter lensing strength is maximized if the angular diameter dis- 
tance to the foreground cluster is about half that of the back- 
ground source galaxies, and to the requirement of large field 
images, in order to cover the full cluster in a single expo- 

, sure. Very few low redshift clusters have been analyzed un- 
til now in detail through weak lensing: Abell 3667 at redshift 
0.055 (Joffre et al. 2000), and clusters observed in the shallow 
Sloan Digital Sky Survey (Sheldon et al. 2001), including Coma 
(Kubo et al. 2007). We present here a weak lensing study of the 
Coma cluster at zcoma = 0.024 based on deep exposures in sev- 
eral bands obtained with the CFH12k and Megaprime/Megacam 
large cameras at the Canada France Hawaii Telescope. The 
J2000 location of the brightest cluster galaxy NGC 4874 is 
RA=12:59:35.67, DEC=+27:57:33.7 . In the following we as- 



* Based on observations obtained with MegaPrime/MegaCam, a joint 
project of CFHT and CEA/DAPNIA, at the Canada-France-Hawaii 
Telescope (CFHT) which is operated by the National Research Council 
(NRC) of Canada, the Institut National des Sciences de I'Univers of the 
Centre National de la Recherche Scientifique (CNRS) of France, and 
the University of Hawaii. This work is also partly based on data prod- 
ucts produced at TERAPIX and the Canadian Astronomy Data Centre 
as part of the Canada-France-Hawaii Telescope Legacy Survey, a col- 
laborative project of NRC and CNRS. 



sume the "concordance model" cosmological background with 
Ho = 70kms"'Mpc"', Q^, = 0.3 and Qa = 0.7, for which 
1" on the sky corresponds to a transverse scale of 0.79 kpc. All 
magnitudes are in the AB system. 



2. Data analysis 

2. 1 . Cluster membership & background sources 

Multiband imaging for this cluster was obtained at CFHT, both 
in the inner parts made of two CFH12k camera pointings in 
BVRI bands resulting in about 0.54 deg^ coverage and with 
Megacam (1x1 deg^ fov) in the u* band (see Adami et al. 2006, 
2007, 2008, for a detailed description of these data). 

We were able to estimate photometric redshifts of all sources 
down to /? < 23 in the area covered by the CFH12k images 
(Adami et al. 2008). This allowed us to separate the popu- 
lation of faint background sources from nearby objects with 
Zphot ^ 0.04 and to eliminate the bulk of the cluster members that 
could bias the weak lensing study (e.g. Medezinski et al. 2007). 
Adami et al. (2008) used the same technique to discriminate be- 
tween z < 0.2 and z > 0.2 galaxies. Here, we push the limit 
to Zphot - 0.04 to increase the number of background sources 
although a more conservative cut at 0.2 would not produce no- 
ticeable changes given the small number of added objects. 

This photometric redshift separation is possible only in the 
innermost 0.54 deg^. However we found that none of our results 
below could be significantly biased due to the extreme preva- 
lence of distant sources in the faint end of the magnitude counts 
(a factor ~ 10 - 20) despite some incompleteness for u* > 24.5 . 
We thus apply a magnitude selection 23.5 < m* < 25 everywhere 
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in the field but we complement the selection by removing the few 
objects in the inner regions that are known to be at Zphot < 0.04. 

A clear advantage of lensing by nearby deflectors is that the 
geometric distance factor Dls/^s {e.g. Bartelmann & Schneider 
2001), which is the ratio of angular distances between deflector 
and source Dls and to the source Ds, is almost constant and 
close to unity for most of the faint background sources. Given 
the low redshift of Coma, even a source at z ~ 0.2 will experi- 
ence lensing effects with Dls/^^s - 0.87. In order to check this 
small redshift dependance, we used the CFHTLS-deep photo- 
metric redshifts (Ilbert et al. 2006) measured in four independent 
1 deg^ fields to infer the redshift distribution of sources having 
23.5 < u* < 25. The mean and median photometric redshifts 
are 0.87 and 0.92 respectively and the sample contains a negli- 
gible amount of low redshift (z < 0.04) galaxies. We estimate 
the mean w - (Dls/^^s) averaged over the whole population 
of sources to be vv ^ 0.962. This is the value we will consider 
below for further weak lensing mass calibrations. 



e,= (l.009±0.026)eu 




Fig. 1. Corrected ellipticities of galaxies measured in Megacam-M* and 
CFH12k-/ bands in the inner parts of the field. A close consistancy is 
found withe; = (1.009±0.026)e„. with a ~ 0.77 correlation coefliicient. 



2.2. Shear measurement 

In Adami et al. (2009), we presented measurements of the sys- 
tematic orientation of faint Coma cluster galaxies. The small 
size of these galaxies required a careful deconvolution of the 
Point Spread Function (PSF) smearing effect which was done 
with the same procedure we use here. It is based on the KSB 
technique (Kaiser et al. 1995) and the pipeline has been applied 
to ground-based imaging from CFH12k and Megacam cameras 
(e.g. Gavazzi et al. 2004; Gavazzi & Soucail 2007). We used the 
u* channel to measure shapes because it is the one with largest 
spatial coverage, and also the most homogeneous since the im- 
ages in the other bands are built from two adjacent pointings. 
The median u* band seeing FWHM is r.'02. Despite this rela- 
tively poor value, we were able to successfully correct for PSF 
smearing. 

The quality of PSF anisotropy corrections can be seen in the 
left panels of Fig. 1 in Adami et al. (2009). The good overall 
image quality is reflected in the raw and corrected ellipticities of 
stars. The low scatter in the corrected ellipticities, 0.0028 rms, is 
also a good assessment of the controled systematics. 



Based on mock images designed to simulate Megacam data 
for the CFHTLS survey and kindly made available to us by 
B. Rowe, we were able to assess the reliability of our shear mea- 
surements, especially for the isotropic part of the PSF smear- 
ing effect. These simulations are similar to STEP2 simulations 
(Massey et al. 2007) as they capture the complex shape of 
sources by the use of shapelet models. Following notations of the 
STEP shear measurement project (Heymans et al. 2006; Massey 
et al. 2007), residual additive calibration errors are found to be 
c < 0.003 and multiplicative m < 0.05, which is significantly 
smaller than the statistical errors present in our data (see below). 

Because redder bands are usually used for weak lensing stud- 
ies, we checked that PSF-corrected ellipticities of central objects 
measured in both Megacam-w* and CFH12k-/ bands are consis- 
tent. This is found to be the case with great accuracy as shown in 
Fig. 1 . We checked that the ellipticity difference e/ -e„. (for both 
components 1 and 2 taken together or considered individually) 
does not correlate with either e„., u* magnitude, / magnitude or 
r/i i,. the M*-band half light radius. 

2.3. Convergence map 

From the source catalogue, we can infer the shear field y(0) and 
deduce the associated convergence field k(6) which is the pro- 
jected surface mass density expressed in units of the critical den- 
sity E„it = c2(47rG)-'Ds/(DLZ)Ls). Here Ds/Dls is replaced by 
the source population average inverse 1 / w which leads to a crit- 
ical density Zcrit - 1-69 x lO'^/zvoMo/kpc^. 
Shear and convergence are related by: 

k(0)^ f K(e-&yy(&)(f&, (1) 

K(6) = ^ (fli-ik)- ^ complex convolution kernel (Kaiser & 
Squires 1993). The shear field is smoothed with a Gaussian 
0(6) oc exp(-^) with 0, = 170" - 80/!7,jkpc. The conver- 
gence field is smoothed by the same filter The convergence map 
presents correlated noise properties (van Waerbeke 2000) char- 
acterised by nbg - 11 arcmin"^, the number density of back- 
ground sources. By bootstrapping the orientation of background 
sources, we estimated this noise level to be cr^ - 0.0134 with lit- 
tle variation over the field (except near edges and masked bright 
stars). 

The top panel of Fig. 2 shows the reconstructed convergence 
map 1x1 deg^ field of view with contours showing a > 5cr 
significance central peak along with various substructures. The 
high physical resolution of ~ 80 h^^ kpc that is made possible 
by the vicinity of Coma somewhat balances the low associated 
lensing efficiency. The position of the brightest cluster galax- 
ies with r < 14 is shown as red crosses. For comparison we 
also show in the lower panel of Fig. 2 the luminosity-weighted 
r-band luminosity distribution coming from bright r < 19 clus- 
ter members that lie on the Red Sequence defined by the rela- 
tion 0.72 <ig-r) + {r- 14) * 0.035 < 0.87. The convergence 
contours are overlaid. For this latter study we used the more ex- 
tended SDSS DR7 photometry, which is better suited for inves- 
tigation of the bright end of the luminosity function'. This panel 
also shows contours of X-ray emissivity in excess of a smooth 
y6-model profile (Neumann et al. 2003). 



' http : //cas . sdss . org/astro/en/tools/ 
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Fig. 2. Top panel: Convergence map for the Coma cluster Green 
contours represent signal-to-noise ratios of 1,2, ...5, corre- 
sponding to K - 0.01, 0.02, . . .0.05. Red crosses represent the 
bright cluster members lying on the red sequence. Bottom 
panel: Gray-scale view of the luminosity distribution of clus- 
ter red sequence members with the overlaid contours in green. 
The blue dotted contours show the excess X-ray emission over a 
smooth yS-model X-ray emissivity map (Neumann et al. 2003). 



3. Radial mass profile 

We now investigate the azimuthally-averaged tangential shear 
profile yriR) which is simply related to the azimuthally-averaged 
projected surface mass density profile by the relation 



^cntYtiR) = AE(-R) = M(< R)/i7rR^) - ^R) , 



(2) 



where we have defined the frequently used rescaled shear AS 
and M(< R) is the projected mass enclosed by radius R. 

We measure the radial shear profile average in circular annuli 
centered on the peak of the convergence map (q'=12:59:39.007, 
5=H-27:57:55.93) which is only 40" east of NGC 4874, i.e. in 
the direction of NGC 4889, which are respectively the brightest 
and second brightest member galaxies. Fig. 3 shows the radial 
shear profile jiiR) out to the edge of the Megacam field of view 
(~1 Mpc). The bottom panel of this figure shows the same pro- 
file once galaxies are rotated by 45°, which is the curl or B-mode 



component of the ellipticity field. In the absence of systematic 
PSF correction residuals, this rotated shear profile should be con- 
sistent with zero at all scales. This is what we observe. 
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Fig. 3. Top panel: Radial shear profile measured in Coma as a 
function of projected distance. The best-fit NEW profile is over- 
laid. Bottom panel: B-mode shear profile showing the negligi- 
ble residuals in the rotated shear component. 



We attempted to fit a radial shear profile as predicted by the 
NFW mass density profile of the form 



p(r) = ps(r/rs) (1 + r/r.T 



(3) 



5.1!4-3 X Wh-l^Ms and caooc 



5.0^11, which 
1.8;^y3Mpc. Fig. 4 also shows the mass- 



coming from cosmological simulations (Navarro et al. 
1997). The corresponding lensing quantities were derived 
by Bartelmann (1996). The two relevant quantities that 
we fit are the mass M200C enclosed in the radius r2ooc in 
which the mean density is 200 times the critical density 
Peril and the concentration parameter C2ooc = '"200c /'"j- This 
implies that - p^nt^c^ / [ln{l + c) - c/(l + c)]. Fig. 4 
shows the constraints we obtain on these two parameters. 
Marginal distributions yield the following constraints^: 
M200C 

corresponds to r20ft 

concentration relation and its Icr dispersion that were recently 
reported by Maccio et al. (2008) assuming WMAP5 cosmo- 
logical parameters (Komatsu et al. 2008). The two-dimensional 
contours are in good agreement with these predictions. Taking 
the conditional p(c2ooc|A^200c) from Maccio et al. (2008) as a 
prior on C2(K)c, we marginalize again over the poorly constrained 
concentration parameter and obtain constraints on the mass 
^200c = 9.7!«;i X lO'^/z-^iMo or again r2ooc = 2.2!0;3Mpc . 

In order to allow a comparison with other mass estimates in 
the literature, we calculate the virial mass Mvii- s Miooc and the 
corresponding concentration Cvii = ciooc, since a density con- 
trast Avir - 100 is better suited for the assumed cosmology. All 
our results regarding both (M200C, C2ooc) and (Mvir, Cvii) as well 
as corresponding r2ooc and rvir values are given in Table 1 

These results are in good agreement with mass estimates in 
the literature. Kubo et al. (2007) performed a weak lensing mass 
estimate of Coma based on the much shallower SDSS data - re- 
sulting in «bg - larcmin"^ - but extending out to ~ 14 /i^^ Mpc. 



' Assuming a flat uniform prior in log M200C betwen 13 and 16 and on 
c'200c between 0.01 and 20 to infer respectively the marginal PDFs p(c) 
and p(M2ooc) 
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Table 1. Summary of mass estimates with and without priors on 
the mass-concentration relation (Maccio et al. 2008). 





no prior 


prior p(c\M) 


M2ooc/10"*/i7^Mo 


5.1+^-^ 


Q 7+fJ.l 


C200C 




3 5+^-^ 

-Q-9 
9 9+0.3 

"^■"^-0.2 


riooc/h-^ Mpc 




M™-/10»A,-JMo 


6 r'-' 

-3(5, 


11 r'^'' 

-6-J 

4 9+1.7 
3 6+'-' 


/i^o' Mpc 





They report a mass M200C = 2.7;^| g x IO'^/i^qMq and con- 
centration index C2ooc = -^.S^Pg' that are statistically consis- 
tent with our estimates. In addition, an earlier X-ray study by 
Hughes (1989) gave Myir = 13 + 2 x h-^^Ms, and is consis- 
tent with our estimate. Mass measurements coming from galaxy 
kinematics are also in good agreement with our mass constraints 
(e.g. The & White 1986; Geller et al. 1999; Lokas & Mamon 
2003). Note however that these latter authors report a concen- 
tration parameter Cyii - 9, slightly larger that our best-fit value 
though statistically consistent. 
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Fig. 4. M200C and C200c parameters constraints of the NFW den- 
sity profile. The M - c relation predicted by simulation and 
WMAP5 cosmology (with its Icr scatter) are shown in red. 



4. Discussion and conclusion 

We have performed a weak lensing analysis of the Coma cluster, 
based on a deep u* exposure and on images in other bands which 
allowed to derive photometric redshifts in order to remove as 
many foreground galaxies as possible from the background pop- 
ulation. We thus obtain a map of the mass distribution in Coma, 
the radial shear profile, as well as the mass and concentration 
parameter. 

The peak of the convergence map is nearly coincident with 
the X-ray center (Neumann et al. 2003) or the brightest galaxy 
NGC 4874. The difference is ~ 27 /i^^' kpc, below the spatial 
resolution of the convergence map. We can also notice the cor- 
respondence between the (a~194.4, 5~27.82) 2cr detection and 
the G12/G14 groups of Adami et al. (2005) and between the 
(a~l94.5, (5~28.03) Icr detection and the northern part of the 
west X-ray substructure (Neumann et al. 2003). Finally, the con- 
vergence map shows a 4cr (q'~195, 5~28) extension overlapping 



the position of NGC 4889, confirming the existence of a mass 
concentration around this galaxy. Other substructures visible in 
X-rays and in the optical (e.g. the NGC 491 1 group) do not show 
up in the convergence map, indicating that we are dealing with 
low mass systems. 

The south-west extension of the convergence map is more 
puzzling. It does not match the light distribution. The groups 
G8 and G9 from Adami et al. (2005) explain part of this exten- 
sion, but are not luminous enough to fully account for this mass 
concentration. Another explanation would be massive low red- 
shift structures (in order to still efficiently induce a weak lensing 
signal and to be spread largely enough over the sky) but that 
would not belong to Coma. The SDSS data used by Gott et al. 
(2005) exhibit such az ~ 0.08 "Great Wall", probably not lumi- 
nous enough though to explain the south west mass concentra- 
tion. Moreover, the magnitude limit of the SDSS spectroscopy is 
not very deep, leaving room for other possible large scale struc- 
tures not appearing in the SDSS sample. This will have to be 
confirmed by a more detailed knowledge of the immediate back- 
ground of the Coma cluster (Adami et al. in preparation). 
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